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Abstract 

We propose a relation between the elliptic SL(N, C) top and Toda systems and obtain a new class of 
integrable systems in a specific limit of the elliptic SL(N, C) top. The relation is based on the Inozemtsev 
limit (IL) and a symplectic map from the elliptic Calogero-Moser system to the elliptic SL(N, C) top. In 
the case when N = 2 we use an explicit form of a symplectic map from the phase space of the elliptic 
Calogero-Moser system to the phase space of the elliptic SL(2, C) top and show that the limiting tops are 
equivalent to the Toda chains. In the case when N > 2 we generalize the above procedure using only the 
limiting behavior of Lax matrices. In a specific limit we also obtain a more general class of systems and 
prove the integrability in the Liouville sense of a certain subclass of these systems. This class is described 
by a classical r-matrix obtained from an elliptic r-matrix. 
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1 Introduction 



We study four integrable systems whose equations of motion have a Lax representation with spectral parameter 
PP, [2J, [3]. We consider periodic and non-periodic Toda chains, the elliptic Calogero-Moser model, and the 
elliptic SL(N, C) top. It was established earlier by several authors that the systems are related to each other. 
In |4J Inozemtsev has proposed a procedure (IL) giving a limit relation between the Toda chains and the elliptic 
Calogero-Moser model. Later, the IL was generalized and used to establish connections between other integrable 
systems. Chernyakov and Zotov have shown in [5] that the IL applied to the SL(N,C) elliptic Euler-Calogero 
model and the elliptic Gaudin model produces new Toda-like systems endowed with additional degrees of 
freedom corresponding to a coadjoint orbit in SL(N,<C). Levin, Olshanetsky, and Zotov [5] have constructed a 
singular symplectic transformation from the elliptic Calogero-Moser system to the elliptic SL(N,C) top. Using 
this transformation Smirnov has shown in [7] that integrable tops on the algebra sl(N, C) are equivalent to 
the TV-particle trigonometric and rational Calogero-Moser systems. The relations between the systems can be 
described by the following diagrams: 

ECM model < — > SL(N, C) top 

|IL (1.1) 
Toda system 

Elliptic CM model i — ► SL(N, C) top 

I I 
Trigonometric/Rational CM model < — > Limiting top 

First goal of this paper is to obtain the Toda chain from the elliptic top. This complements the diagram 
(jl.ip in the following way: 

ECM model < — > SL(N, C) top 

| IL | IL 

Toda system < — > Limiting top 

In order to obtain a new relation we will use a procedure similar to the Inozemtsev limit. The Inozemtsev 
limit is a combination of the trigonometric limit, infinite shifts of particle coordinates, and rescalings of the 
coupling constants. To obtain a limiting system equivalent to the Toda chain it is necessary to combine the 
Inozemtsev limit and the infinite shift of the spectral parameter. Since the spectral parameter of the elliptic 
SL(N,C) top is given on a complex torus T 2 with moduli r, under the trigonometric limit Im(r) — > +oo we 
obtain systems with spectral parameter given on an infinite complex cylinder C/Z. 

In the case of the elliptic SL(2, C) top it is convenient to use an explicit form (|1.13|) of a symplectic map from 
the phase space of the elliptic Calogero-Moser system to the phase space of the top (Subsection II. 3p . Then the 
Inozemtsev shifts of the elliptic Calogero-Moser system coordinates induce the rescalings of the elliptic SL(2, C) 
top coordinates. The equivalence between the limiting systems and the Toda chains is due to the bosonizations 
formulas which follow from the limit of (|1.13l) . 

To derive an explicit form of the map between the phase spaces of the Calogero-Moser system and the 
elliptic top in the case of the elliptic SL(N > 2, C) top is not as simple as when N = 2. That is why we use 
the scalings of coordinates induced by the limiting behavior of Lax matrices and thus generalize N — 2 case. 
Also, the scalings of coordinates satisfy an important requirement, that is the limit of the Poisson algebra of the 
elliptic SL(N,<C) top must define a Poisson structure on the phase space of the limiting system. This Poisson 
structure along with the values of the Casimir functions define the symplectic submanifold, for which there is 
the symplectic map to the phase space of the Toda chain. Equations of motion of the limiting system have Lax 
representation and are equivalent to the equations of motion of the Toda chain. 

Second goal of this paper is to obtain in the limit a more general class of systems and prove the integrability 
in the Liouville sense of a certain subclass of these systems. This class appears under specific conditions on the 
parameters of the limit and contains Toda chains as a special case. It is possible that further study of this class 
will lead to establishing a connection between integrable systems mentioned above and gauge theories. Such an 
approach was developed earlier. For example, in [8], [9] Toda-like systems corresponding to the multi-component 
magnets were studied in the context of the low-energy effective N = 2 SUSY gauge theories. 

Now we will review general facts and notation about the integrable systems under consideration. 
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1.1 Elliptic SL(N, C) top 

The elliptic SL(N,C) top is an example of Euler-Arnold top |10j . The elliptic SL(N,C) top is defined on a 
coadjoint orbit of the group SL(N, C): 

K Tot ={S €sl(N,C), S = a- 1 S (0) a}, (1.2) 

where a € SL(N,C) is defined up to the left multiplication on the stationary subg roup G of S (0) . The phase 
space lZ mt is equipped with the Kirillov-Kostant symplectic form 

uj Iot = Tr (S^daa- 1 A daa" 1 

The Hamiltonian is defined as 

H rot = ~TrSJ(S). (1.3) 
Here we consider a special form of a linear operator J that provides the integrability of the system 

' m + tit 



m,n G {0, ...,N- 1}, m 2 + n 2 ^0, 

where -E^C^j r ) is the second Eisenstein function (see [ 1 1 J ) defined on the complex torus T 2 : C/ (2wiZ + 2^2^) 
with wi = 5, t = u)2/u)\, and s mn are coordinates in the sin-algebra basis T mn (see[A|). 
The equations of motion can be written in the Lax form |12| : 

J T rot 

=N[L Iot ,M Iot ] . (1.4) 
Factor iV in (| 1 .4[) comes from the definition of Lax matrices in the sin-algebra basis from [3] 



L«* = £ SfnnV , [™] (z)T mn , p [™] (z) = e (-^) 



m,n 

nz 



m + nr 

-,z 



N 



[1.5) 



where e (z) = exp (2iriz) , i = \J — 1, and <ji> is a combination of theta-functions (see[B]). The Lax matrix satisfies 
the properties of quasi-periodicity: 

U ot (z + 1) = T 10 L rot (z) Tfo 1 , L rot (z + r) = T iL rot (z) T^ 1 . (1.6) 

Consequently, Tr (Z/ ot {z)) k are doubly periodic functions with the poles of order up to k, and thus they can 
be expanded in the basis consisting of the second Eisenstein function and its derivatives: 

Tr (L lot {z)f = H kfi + E 2 (z) H k<2 + E' 2 (z) H k>3 + --- + E^ 2) (z) H k>k . 
In this way we obtain the Hamiltonian (|1.3p 

H rot = ^H 2 ,o = ^Tr(L mt f - ±TrS 2 E 2 (z, r). (1.7) 



Poisson brackets for variables s mn are defined by the commutator [T a b,T c d] (|A.1[) of basis elements T a b and 

T c d 

rr 

8 a +c,b+d- (1- 



r ft 

{sab, s cd } = 2i sin — (6c - ad) 



Then transition to the standard basis (|A.2|) gives us 

Ski} = N(S k jSu - SuSkj)- (1.9) 

Linear brackets (|1.8[) , (|1.9p can be written in terms of the Bclavin-Drinfeld classical elliptic r-matrix r(z) 
US, H3, IE]- Namely, 

{L?* . i 2 ot (%)} = [r (*i - *&) , L{ ot ( Zl ) + L 2 ot (z 2 )] , (1.10) 



3 



where 



Li (z) = L (z) (g) Id, L 2 {z) = Id® L (z) . 
The classical r-matrix is defined by 



r (z) = -y^tp\ (z)T mn ®T- r 
^— ' L n . 



(1.11) 



Equation (|1.10|l implies the involutivity of the independent coefficients Hk, z - Therefore, there are N(N + 
l)/2 — 1 independent integrals of motion. Note that Hk,k, k G {2, . . . , N}, are the Casimirs corresponding to 
the coadjoint orbit (|1.2I) . 



1.2 Elliptic Calogero-Moser system 

The elliptic Calogero-Moser system (CM) was first introduced in quantum version [16J, [17J. The elliptic CM 
system is defined on the phase space as follows 



n 



with the canonical symplectic form 



CM 



{N N ~j 

(u,v), ^^ = 0, ^2^ = 0} 
i=l i=l J 



uj^ m = (dv A du) . 



The corresponding Hamiltonian is defined via 

N 2 
i— 1 i>j 

The equations of motion denned by the Hamiltonian have the Lax representation 

dL CM 



dt 



= [L CM ,M CM ] 



where the Lax pair can be chosen in the holomorphic form in order to construct the connection between the 
Calogero-Moser system and the elliptic SL(N, C) top [6j 

Lf 3 M = SijVi + m (1 - Sij) 4> (ui - Uj, z) , 



M 



CM 



d<p(u, z) 



dv 



1.3 Connection between the Calogero-Moser system and the elliptic SL(N,C) top 
The connection mentioned in the headline was established in [6 J in the form of a singular gauge transformation 

L*°\z) = Z(z)L CM (z)Z- 1 (z). 

This transformation leads to the symplectic map 

^CM^^rtrt ( U ,v)l-»5. 

In the case when N = 2 this map has the form 



(1.12) 



0oi(O)M2«) 
Sal =~ V *'(0W2«) ~ m \ 

lo (O)6 lo (2u) 
SlO = v — — — ; h m- 



tf'(0)tf(2u) 



0^(0) 9 00 (2u)8 w (2u) 
»oo(0)eio(0) $ 2 (2u) : 

gyp) 6 oo (2u)0 ol (2u) 
>oo (O)0oi (0) $ 2 (2u) ' 



(1.13) 



_. gpo (0)g O o (2m) . ggn(O) 0w(2u)6 01 (2u) 
SU W ^(0)i?(2«) 1W lo (O)0 O i(O) ^ 2 (2 U ) ' 
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1.4 Toda systems 

Periodic and nonperiodic Toda systems of N interacting particles in the center of a mass frame are defined on 
the phase space 



{N N ~| 

(U,V), 5>,;=0, £>i = 0L 
i=\ i=l ) 



with the canonical symplectic form 

uj t = (dv A du) . 



The Hamiltonian of the nonperiodic system is 

N A — I 

2 



H AT = \ J2 vj + 4tt 2 M 2 £ e(u l+1 - «,), 
' »=i i=i 

and of the periodic system has the form 

iV JV 

»=1 i=l 



Ml. 



The equations of motion for both nonperiodic and periodic Toda systems can be written in the Lax form Q18J. 

^L AT = [i AT , M AT ] , A L PT = [ L PT j M PT] _ 

One can obtain Toda systems by applying the Inozemtsev limit to the Calogero-Moser system [4|. 

2 Elliptic SL(2, C) top via the Inozemtsev limit 

The main idea of the technique under consideration is to treat elliptic SL(2, C) top coordinates as functions of 
coordinates (u, v) of the elliptic Calogero-Moser model and apply Inozemtsev shift of (u, v). 

2.1 Periodic Toda system from the elliptic top 
2.1.1 Limit of the Lax matrix and Poisson algebra 

To obtain the periodic Toda system we combine the shift of coordinates u = U + r/4, scaling of the coupling 
constant M = mq^ (q = e(r)), the shift of the spectral parameter z = z + r/2, and the trigonometric limit 
q — > 0. Then from (|1.13[) we derive 

iv _ , i . 

sio = \-0(q*), 

Mcos(2ttJ7) wsin(27r/7) n , i, 

*oi = \ + 0{q*), 

q-i it 

Msm&irU) vcos(2-nU) i, 
ail = 1 + O(g'). 

qi 7r 

Coordinates of the limiting top are scaled coordinates of the elliptic SL(2, C) top 

~ iv 

s w = lim s w = , (2.1a) 

q-*0 7T 

sbi = lim s i9^ = M cos(27rf7), (2.1b) 
sn = lim sng3 = -M sin(27r?7) . (2.1c) 

q— >0 

Scaled coordinates (|2.1ap - ()2.1cll form an algebra which arises via contraction of si (2, C) algebra 

{siojSii} = 2is i, {sii,s i} = 0, {s i,sio} = 2is"u. (2.2) 
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Also, formulas (|2.1ap - (|2.1cp define a symplectic map from canonical coordinates (U, v) to the coordinates 
of the limiting top. Such formulas are known as bosonization formulas. 
The following condition defines the symplectic leaf: 



~2 i~2 

s oi + s ll 



const = M 2 , 

and originates from the Casimir function of the elliptic SL(2, C) top 



'01 "T" l, 10 



Sin + S\i — const — to 2 . 



Taking into account the behavior of the function tp [™] (z) (IB.8|) (see|B|) we can write down the limiting Lax 
matrix 



L Iot = Air 



sqi sin(7rz) — sn cos(7rz) \ 



y Sqi sin(7rz) + sn cos(irz) 



2.1.2 Limiting equations of motion and bosonization 

We will use formula (|1.3p for computing the limit of Hamiltonian 



H rot = _ 



where 



J10 = lim J10 = 7r 

q-S-0 

J01 = lim J a iq~^ = - 

q->0 



in 2 , 



Jn = lim Jn<7 



8tt 2 



The series expansion of / [™] (z) ()B.9[) (see|B]) leads to the second Lax matrix 



4(s i +isu)e 
-sio 

Equations of motion (|1.4p preserve the same form in the limit: 

dL rot 



4(s i -isn)e ( -- 



dt 



l^rot^rotj 



L rot )M rot 



Using bosonization formulas (|2.1ap . (|2. lb[) . and (|2. lc|) one can obtain the periodic Toda system 

H rot -> H PT = v 2 + 8M 2 ir 2 cos (4ttZ7) , 

v AttM sin (n (2U + z)) 

-47rMsin(7r(2C7-J)) -v 

\ 



V 



PT 



M rot ^ M PT = 



-mv 



4:n z Me -U 



An 2 Me U - - 



J 
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2.2 Nonperiodic Toda system from the elliptic top 
2.2.1 Limit of the Lax matrix and Poisson algebra 

Here we will use another shift of coordinates u = U + t/8, another scaling of the coupling constant M = mq 1 / 8 . 
the same shift of the spectral parameter z = z + r/2 and apply the trigonometric limit q — >• 0. That gives us 
the following formulas for the elliptic SL(2,C) top coordinates 



soi 



su 



sio = +0((?«), 

Me(t7) ive(U) 
' i ' i 
2g* 27rqs 

iMe(J7) ue (f7) 



2qi 



2nqt 



O(l), 
O(l), 



and for the coordinates of the limiting top 



sio = hm sio = , 

9^0 7T 



soi = lim soig 4 = -Me(Z7), 
su = lim sn<7^ = -Me(Z7). 

g->0 2 

These coordinates form the same algebra (|2.2[) as in the case of the periodic Toda system. 
The limiting Lax matrix has the form 



47T 



tSw 



soi s 



i sin(7rz) — su cos(7rz) \ 



y soi sin(7rF) + su cos(7r?) 



rSio 



In this case the symplectic leaf of the elliptic SL(2, C) top turns into 

0. 



6 1 T S X j 



(2.3) 



2.2.2 Limiting equations of motion and bosonization 

Since we have the same shift of the spectral parameter and the same scaling of coordinates s mn as in the periodic 
case, J mn have equivalent limits 

Jio = TT 2 , Joi = — Ju = — 87r 2 . 
Also, limiting Hamiltonian and the second Lax matrix acquire the same forms 



h = - r 01 j 01 + s{ j w + 7^ Jn , 



M xot = Tl Z 



SlO 



4(s i +isn)e 
-sio 

The limiting Hamiltonian can be simplified on the symplectic leaf (|2.3p as follows 



4(s i -isn)e ( -- 



2s i Jqi + s w Ji< 



The equations of motion have the Lax representation 

dZ rot 



dt 



{jj rot ,Z rot } 



L rot M rot 
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Bosonization formulas transform the limiting top into the nonperiodic Toda system 



L rot ^ ^AT 



v 2 + 4M 2 n 2 e(2U), 

v -2mMe(U- 
2mMe(U - f ) 



M rot -> M AT = 



- r 




47r 2 Me(C/-f) iTTt! 



3 Elliptic SL(N > 2,C) top via the Inozemtsev limit 

In this section we consider a limit that is a combination of the shift of the spectral parameter z = z + t/2, 
the scalings of coordinates, and the trigonometric limit Im{r) — > +00. The searings of coordinates are defined 
by the limiting behavior of the Lax matrix and are not derived from the symplectic map (|1.12[) as in the case 
N = 2. Also, the scalings of coordinates satisfy an important requirement, that is the limit of the Poisson 
algebra of the elliptic SL(N, C) top must define a Poisson structure on the phase space of the limiting system. 
This Poisson structure along with the values of Casimir functions define the symplectic submanifold for which 
there is the symplectic map to the phase space of the Toda chain. Equations of motion of the limiting system 
have Lax representation and are equivalent to the equations of motion of the Toda chain. 

3.1 Periodic Toda system from the elliptic top 
3.1.1 Limit of the Lax matrix and Poisson algebra 

In order to determine the exact scaling of coordinates we need to expand the function <p (z) as series in q, 
where q — e 2mT (see [Bj . We obtain 















. n - 







-™{2N) Sm hv) + ° (1) ' 



TIZ Tfl 



0(qi) 



n = 0, 



N 

0<n<—, 



N 



(3.1) 



2 ' 



,'N-n. 
— TV — Z ' ^ 2 J ' V 



oyq' 2 2JV J, — < n < 



N. 



Since the Lax matrix for the periodic Toda system can be written in tridiagonal form, the following substi- 
tution is reasonable 



1 - 8(n) 



e {0,...,N- 1}, ra 2 + «Vo, 



g{n) = 

This gives us the limiting matrix U 



2N 



N-l 



(3.2) 



\2NJ 



1 / m\ 
\~N / Sm0 °~ 



Af-l r 



m=0 



z m \ ~ m I z \ 

e ( — — + — I s rn \T m i — e I — ) s TO ,jv-i-im,j\ 1 



(3.3) 



where 



1 n = mod TV, 
n ^ mod TV. 

As one can see, coordinates s mn , 1 < n < A — 1 , are not present in the adduced matrix and hence in the 
Hamiltonian. We will show later that Hamilton equations for these variables can be integrated in spite of the 
fact that their dynamics are separated from the Lax representation. 
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After scaling (|3.2p . we obtain the contraction of Poisson algebra (jl.8p in the limit q — > 

{l ab , I cd } = 2i sin (6c - ad)] q^ +9 ^-^ b+d ^ a+c , b+d . (3.4) 

Therefore, scaled coordinates s mn with the Poisson brackets form an algebra in the limit of q — > provided that 

Vfc,n: g(k) + g{n) - g(k + n) ^ 0. (3.5) 

If g(n) = (l — S (n)j / (2N), then (|3 . 5p is trivial and we can write down all nonzero brackets corresponding 
to the equality in (|3.5p 

{sao,s c d} = -2isin ^adj s a + c ,d- 

It is convenient to use the standard basis further. In this basis substitution (|3.2[) and the Lax matrix turn 
into 

Si^Syq-tM, &(i,j) = ^ i , i,je{l,...,N}, 



+27rLS i+1 , i+2 e S (j - i - 1) - 27rLS M -ie l-j 5 (j - 1 + 1) . (3.6) 

From now on, indexes of S belong to {1, . . . , N} and satisfy the properties of periodicity 

SN+i,j = Si : N+j = Sij i,j € Z. 
From (|1.9[) we obtain the following nonzero Poisson brackets for coordinates in the standard basis: 

{Su, Sjk} = N(SjiSik - SikSij). (3.7) 
Now, it can be easily seen that Casimir functions are 

N 

Si lt i 2 Si 2< i 3 ■ . . Si k . h Vj ^ I ii^ij, 2 ^ k ^ N. (3.8) 
In (|3.8p there are iV + 2 independent functions depending only on variables that form the Lax matrix ()3.6j) 

N N 

Su, JJS'j^-i, Si^+iSi+i^i i € {1, . . . , N}, 

i=l i=l 

Also, in (|3.8p there are N (N — 3) Casimir functions independent as a functions of variables which are not 
included in the Lax matrix 

Thus, on the symplectic submanifold with nonzero values of quadratic Casimir functions S^+iSi+i^ variables 
which are not included in the Lax matrix are the following functions of variables included in the Lax matrix 

k 

Si+k,i = const JJ Si + j,i+j-i 1 < t < iV, 2<fc<A r -2. 
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3.1.2 Limiting equations of motion and bosonization 

The limit of Hamiltonian, as it follows from (| 1 . T[) and the fact that TiS 2 E 2 (z, r) — > 0, depends only on variables 
contributed in the Lax matrix 



H 1 



m,n— 1 k— 1 



k(n — m) 



N 



+47T 2 



Using the series expansion of / [™] (z + r/2) fsee lB]) 

' -^ 2 sin- 2 (tt^) +o(1), 



" m ~ 










"I) 


. n - 







= < 



(£)-(- 
( 



g2JV + o ( g2« 



iV 



<• i —-z ) qs +o[q> 



, , m „\ / nz\ 3 

-47T e h z e , 

N J \ N 



one can obtain the second Lax matrix 



k 
N 



n = 0, 



< n < 



3N 



3N 



4 ' 



) e (~^)' ? ^ 1 ™' ) + ('7^ 1 "'')' < n < ^ 



(3.9) 



2V-1 TV 



m=lfe=l ^ ' 

and ensure that the equations of motion can be written in the Lax form 

d 



z 

Nj ' 



L xot M rot 



dt 

Those variables that are not included in the Lax matrix (|3.6[) have simple dynamics 



(3.10) 



(3.11) 



-^sv, = 47r 2 sv, e E s ' mm sin ( 77 

TV 

k 



N N-l 



x 1 — cos 2ir 



1)1 = 1 A:=l 
-1 



k(j — i)\ ( k(i + ? — 2m) 

1 sin tt— 4^ I x 



iV 



N J \ N 
1 < (j - i) mod N < N — 1, 



and with other coordinates allow bosonization formulas 



(3.12) 



5 N < \ 

Z7T1 

Si ti+1 =MNe( Ui ), 

S l+hl = MNe(-Ui), 

/i+k-1 \ 

Sj,»+Ai = Ci )i+ fee I E u " ) 2 ^ fc < TV - 2, c^i+fc = co?isi, 

\ n=i / 

where u, v are canonical coordinates 

{vi,Uj} = Sij i,j € {l,...,iV}, 



and 



E u i = Qj E v * = °' 



(3.13) 
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Let us show that for u, v we have dynamics of the periodic Toda system in the center of a mass frame. 
Substituting (|3.13p into Hamiltonian (|3.9|) and Lax matrices (|3 . 6[) . (|3.10|) we obtain 



N 



H rot = N 2^2^+ 4n 2 M 2 N 2 £ e (u i+1 - u % ) = N 2 H 

i=l i=l 

where H PT has the form of periodic Toda Hamiltonian, 



PT 



L rot = 2mMN 



2?riM „ 



V 2 



2iriM 



-e(-ua + tt) 







«3 



AT 2th Af 










vn-i 

2-k\M „ 
-e(-«jv-i + ^) 



-e(-u N + — ) » 





7V y 



27riAf 



Af, r 



. AT-l N 



2 / m\ / m(i — fc) \ 
sin 1 7T— ] e | — | (v fe _i - Ufc) dij+ 



N 



N 



+4n 2 MNe u l+l - — )S (j - i - 1 



N 



After the gauge transformation we have 



(3.14) 



iV' 



^ e ( ^ ) n e ' 



fe=i 



and Lax matrices take the standard form 

2wiM 
-1 



L Iot = 2iriMN 







e(ui — mat — z) 



e(u 2 - ui) 

2?riM 
-1 








e(w 3 - u 2 ) 

2tt\M 










vn-i 
2it\M 

-1 



-e(2) \ 



e(wAT - ujv-i) 



27riM 



(3.15) 



47r 2 MiV 



\ e(ux -u N -z) 



e(u 2 - wi) 










e(it 3 - u 2 ) 







e(uN - un-i) 
/ 



3.2 Nonperiodic Toda system from the elliptic top 
3.2.1 Limit of Lax matrices and Poisson algebra 

To obtain the Lax matrix of the nonperiodic Toda chain we are going to consider another substitution 



c. . a -a(i,j) 



1 - St. 



— SuSjN 

2 



2N 



i,je{i,...,N}. 



(3.16) 
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Hence, the contraction of Poisson algebra p. 91) in the limit q — > takes the form 

{S lJ: S kl } = Nq'M+'M (tuSkjq-'W - S kj S u q- B{u) ) , 
and the scaled coordinates form an algebra in the limit q — > if 

V«,j, k Q(i,j)+Q{k,i) - &(k,j) ^ 0. 



As it can be easily seen the above inequality is valid for defined in (|3.16|) . Upon the limit g->0we have 

nonzero brackets Q3.7[) and the following Lax matrices: 



„ . N N-l 



N 

m=l fc=l 

-27riS , j + i, J ; +2 e 



k(i — m) 
N 



k 
'N 



- 1 



j^J 5(j-i- 1) - 27riSi,i_ie ( — ) 5 lJ+1 , 



(3.17) 



y N-l N 

r 

— < 

N 



m\ / m(i — fe) 



m=l fc=i 



AT 



+47r 2 £(.?-i-l)S l+M+2 e 




iV 



As long as the algebra has the same limit as in the periodic case we have Casimir functions 
there are N + 1 independent functions formed only by variables contributed in the Lax matrix 



N 



N 



/] Sjii J"J S'i,t+l'Si+l,i i G {1, . . • , iV — 1} 



i=l 



(3.18) 
But now 



and N(N — 3) + 1 Casimir functions independent as a functions of variables which are not included in the Lax 
matrix 



Si,nSn,i, 

k 



5, 



1 < i < AT, 2<fc^A^-2. 



3.2.2 Limiting equations of motion and bosonization 

The limit of Hamiltonian after substitution (J3TT5J) is 



I ^ 2 N N ~ 1 

-Tr(L rot ) 2 = — — ^ ^ S mm S nn e 



k(n — to) 
TV 



AT 



-47T 2 



1 — cos 2tt 



N 



+ 



(3.19) 



The equations of motion can be written in the Lax form 

— Z rot = {H xot ,L Iot \ = N 
at 



L rot ;M rot 



These equations imply simple dynamics (13. 12)) for variables that are not included in the Lax matrix. And for 
all coordinates of the limiting system there are bosonization formulas 



Su = -Vi), ie{l,...,N}, 

ZTTl 



H,i+1 



= MNe{ui), ie{l,...,N}, 
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5 i+ i,i = MNe(-Ui), ie{l,...,N-l], 
Si,n = const e (— ujv) , 



n+k-i 



S iyi+ k = Ois + ke I u n\, 2 < k ^ N - 2, c M+fc = const. 



(3.20) 



Canonical coordinates u, v have dynamics of the nonperiodic Toda chain in the center of a mass frame. 
After substituting (|3.20[) into Hamiltonian (I3.19P we obtain 

N 2 N- 1 

jjrot ^ N 2J2^+ 4n 2 M 2 N 2 e ("if l - «i) = N 2 H AT , 



where H AT has the form of nonperiodic Toda Hamiltonian. Lax matrices take the usual form under gauge 
transformation (|3 . 1 5|) mentioned in the periodic case 



/ 



Vl 



L rot = 2mMN 



2-kiM 
-1 







e(ui -Un-z) 



e(u 2 - Ui) 

2tt\M 
-1 








e(u 3 - u 2 ) 

2niM 










VN-l 

2ir\M 
-1 







e(u N ~ u N -i) 



2niM 



M 1 



A-k 2 MN 



\ e(ui 









- UN 



2) 



e(u 2 - Ui) 










e(u 3 - u 2 ) 







e(u N - u N -i) 




3.3 More general class of limiting systems 

In the previous Subsections we have considered substitutions of variables (|3.2p and (|3.16p , which after applying 
the Inozemtsev limit lead to the Toda chains. It turns out that the substitutions mentioned above are not the 
only possibility to provide the integrable systems in the limit. We will consider the following generalization: 



m,ne {0,...,N- 1}, m 2 + n 2 ^0, 

k N 

, < k < p < — , 

2N' y 2' 



P 



l{i)= { P<k<N-p, 



(3.21) 



N -k 
2N ' 



N~p^k<N, 



where i G Z, k = i mod N, and prove the integrability of the limiting systems in the case when N and p are 
relatively prime. 



3.3.1 Limit of Lax matrices and Poisson algebra 

Scaled coordinates with Poisson brackets (|3.4j) form a Poisson algebra in the limit q — >• provided that the 
following condition is valid 

\/k,n: g(k)+g(n)-g(k + n)^0. (3.22) 

For g(n) under consideration (|3.22[) is proved in Id 
Nonzero limiting brackets have the form 
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{s a o, s c d} = -2isin (-^ arf ) s a+c4 , 

{sab, Scd} = 2isin ^ (be - ad)j s a+c . b+d 
(0 < a < p) A (0 < b ^ p) A (0 < a + b < p), 

{sab, s cd } = 2isin ^ (6c - ad)j s a+c , b+d 

{N-p^a<N)A(N-p^b<N)A{2N-p^a + b< 2N), (3.23) 

or in the standard basis 

{Su,Sjk} = N(Sji5ik — SikSij), 
{Sij, Ski} = N(SkjSu — SuSkj) 

(0 < (j - i) mod N ^ p) A (0 < (I - k) mod N ^ p) A (0 < (j + I - i - k) mod N < p), 
{Sij, Ski} — N(Skj5u ~ SuSkj) 

(N-p^(j- i) mod N < N) A (N - p < (I - k) mod N < N)A 

A(N~p^ (j + l-i-k) mod N < N), (3.24) 

where as usual A stands for "and". 

Formulas (|3.23|) (or (|3.24D ) imply that the limiting Poisson algebra is solvable. Thus, there is no general 
method to construct all Casimir functions, but in the special case when N and p are relatively prime we able 
to present the whole set of independent Casimir functions. 

At first we are interested in Casimir functions in general case when p < N/2. Since elements Si^+k, k G 
{p, . . . , N — p}, have nonzero brackets only with coordinates Su, Si+k,i+k, we obtain the second and iV'th order 
Casimir functions 



S iti+ kSi+k,i kelp,..., y j, i€{l,...,N}, 

N 

II^M+fc k£{p,...,N-p}, ie{l,...,N}, (3.25) 

i=l 

where [x\ is the floor function of x. 

In the case when N and p are relatively prime we can also construct the following N(N — 2p~l) independent 
Casimir functions 

J\ Si+U-l)p,i+3P ) $i+kp,i P< kp mod N < N - p, k,i £ {1,..., N}. 

If p < (N — 1) /2 it is convenient to consider two disjoint subalgebras. The first one is generated by the 
variables 



Si = {Sij, p< (j - i) mod N < N - p\ 



with simple dynamics, which we are going to show further. Equations of motion for the elements of the second 
subalgebra 



So 



{Sij, (0 < (j - i) mod N < p) or (N - p < (j - i) mod N < AT) j 



have Lax representation, so we need to obtain the number of independent Casimir functions from universal en- 
veloping algebra of this subalgebra. Since Casimir functions lower the dimension of the symplectic submanifold, 
we need to majorize the degeneration factor of Poisson tensor nWKkl) (S 2 ): 

{F (S 2 ) , G (S 2 )} = (S 2 ) d {ij) F d [H) G, d m = -§-. 

oSij 
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Formula (|3.24|) implies that Poisson tensor can be represented as a block lower matrix with respect to 
antidiagonal (|D~Tj) . Rank R of this (2p + l)N x (2p + l)iV-matrix satisfies the condition R ^ 2p(N - 1) (|DT2|) 



N 



(see [Dj , which restricts the number of the independent Casimir functions up to N + 2p (here we treat Su 

i=l 

as a Casimir). 

Linear brackets (|3.23p and (|3.24|) can be written in terms of an r-matrix. Namely, 

[IT (Si) , IT (z 2 )} = [r {z x - z 2 ) , L[ ot (zi) + IT (z 2 )] , (3.26) 

where 

Li (z) = L (z) (g) Id, L 2 {z) = Id®L{z) 1 
and matrix r (z) is the limit of elliptic r-matrix (II. lip 



N-l 

r(zi-Z2)= lim r (z x - z 2 ) = vr V] (cot — - cot (tt(zi - z 2 )) ) T m0 <g> T_ m . - 

7m(r)— > + oo z — ' V iV / 



m— 1 

N-l / ,~ ~ - x N-l 



-* sin- («(zi - z 2) ) e ( £ e ( ) £ T mn ^ T. 



2 / ^ V TV 

7 n=l m— 



m,— n- 



Explicit expression for the elements of r (z) 



N-l / /• • \\ 

MO.WiOt 2 ) = n <L e I m J (, cot -jv" + 1 )°U s iih + 

m=l ^ 

2?riiVe(z) / - i) mod iV x 



l-e(2) V ^ 



where we exclude one summand proportional to Id ® Id. 

Substitution (|3.21[) preserves in the limit coordinates of the Lax matrix with respect to the following sin- 
algebra basis elements: 



T mn , n £ {-p, . . . ,p}, 



which gives 



N-l 



L rot = n £ e ^ gin -l ^ ~ moTmQ 
m— 1 

p N-l r , ,~ v 

Ex -> / iz m \ „ 

; — i ™ — n L \ / 



(3.27) 



1=1 m=0 

In the standard basis the Lax matrix acquires the form 



= 2vri^ ( S i+ i, i+1+l e (--£) S (j - i - I) - S w e )S(j-i+ I) 



i=i 



It is convenient to use both the above and gauge transformed form of the Lax matrix 

iz 



= 9 n ot g-\ 9ij = %e ( ^ ) . (3.28) 



Denoting w — e (z) we obtain 
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/ II 523 



-S21 h S34 



Lf = 2ni 



52,2+p 



-wSi. 



N-p+l 



S3,3+p 



— S32 



-wS. 



pN 



-s, 



5p+2,p+3- 



p+2,2 







->Sp+2,p+l lp+2 



... 

Sn-p+2,2 











w— p 



p+l,iV-p+2 



p+1,1 



-5*jv- 



p+l,JV— p 



—Sn,n- p 



In-i Sn.i 
-Sn,n-i In 



where 



N N-l 



hi Lu - w E E s 



„,,„<■ I feil— ?Zij ) [ e I 1-1 



m=l fe=l 



TV 



N 



and consequently 



^■ji — li— 1 



(3.29) 



3.3.2 Lax representation of limiting equations of motion 

Limiting Hamiltonian and the second Lax matrix have similar structures as in the case when the limiting system 
is equivalent to the periodic Toda chain (see Subsection 13. ip . namely 



2 V J N 

p N 

-Ait 2 ^ ^ Si + i.iSi + i.i + i + i 
1=1 »=i 



H«* = lTr(Z™*Y =~ Y^S mm ~S nn e ( MlL^l ] [ , ( , )S j >, 



m.n—1 k—1 



N 



(3.30) 



2 N N-l 



m=l fc=l 



fc \ / k(i — ni) 



N 



- 47r2 E^+i^+i e (-^) 



S(j -i-l). 



When p = 1 these formulae turn into (|3.9[) and (|3.14[) , respectively. 
After the limit equations of motion also have Lax representation 



dt 



L rot M rot 



(3.31) 
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In the case when p < (N — 1) /2 there are variables Si which are not included in the Lax pair. Hamilton 
equations for these variables are 



— Sij = 4n 2 S t j E E ^" 

71= 

k 



N N-l 



k(j — i)\ ( k(i + j — 2m) 
sin I 7r v — ' sin it— — I x 



x 1 — cos 2ir 



m=l fe=l 
-1 



N 



N J \ N 
P < (j - i) m °d N < N - p. 



On the symplectic submanifold with nonzero values of the second order Casimir functions Si } i+ p Si+ P} i, i G 
{1, . . . , N} in the case when N and p are relatively prime variables Si are the following functions of variables 
S 2 

k 

= const Y[ Si+jp,i+(j-i)p, P < kp mod N < N p, k, i € {1, . . . , N}. 



Thus, if we solve the equations of motion (|3.31[) . we immediately obtain the solutions of the equations of motion 
for variables Si. 

3.3.3 Integrability 

The Lax operator of the elliptic SL (N, C) top satisfies properties of quasi-periodicity (|1.6[) . Namely, 

U ot (z + 1) = T 10 L rot (z) T w \ U ot (z + t) = T m U ot (z) T^ 1 . 
After taking the trigonometric limit Im (r) — > +oo the Lax operator has only one quasi-period 

Z rot (z + l)=T 10 Z rot (z)T w \ 

Since Tr (^L rot (z)j are periodic functions in z, they can be expanded in Fourier basis {e (jz) = j G Z}. 

From the gauge transformed Lax matrix U° it follows that there are finite number of nonzero coefficients in 
this expansion. 



Proposition 3.1. The trace of the k-th power of the Lax matrix has the form 

k AI 

2r(i rot (z)) = H k 3 w j , where M = 



j=-M 



kp 
iV 



w — e(z) . 



(3.32) 



Proof. Replacing e (z) by w in formula (|3.27p we obtain a convenient form of the Lax matrix 



N-l 



p N-l 



^ = - E « ^T m0 +2^iEE [e (£) 



^ N s 7nn T 7nn w N s m n T rn n 



m— 1 

p N-l 



n—1 rn—0 



III Is 

n——p m=0 



Then 



(V ot (T)) = y ■■■ E w J2 ' =i " Yl c ^ m ^ n ^ Sm ^ Tm ^- 

mi,ni mfc,nfc i— 1 

By the properties of T mn (see [AJ the following condition holds 



(3.33) 



As rii € {— p, . . . ,p} for any i, we derive the second condition 
(13321 . 



E< 



(3.34) 



^ kp, which along with Q3.34[) implies 



□ 
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Representation (|3.26|) of linear brackets (|3.23p and (|3.24p provides us with the involutivity of the coefficients 
H kj . 

Proposition 3.2. \21\ . The coefficients H^j are in involution, i.e., 

{Hk in , Hk 2 j 2 } = 0- (3.35) 
Proof. Exactly as in the case of the elliptic top we have 

Tr (Z rot (zi)) fcl ,Tr (Z rot (z 2 )) ^ \ = Tr j (Z rot (?i)) * , (i™*^))*" 

Then, it follows from (|3.26p that these functions Poisson commute. Using the expansion (|3.32[) we get the 
involutivity of the coefficients (|3.35p . 

□ 

In particular, all functions H^j Poisson commute with the Hamiltonian (|3 . 301) . Moreover, among coefficients 

Hf^j we have Casimir functions of the form H k ^ ±j, k(j) = j^- . j g {1, . . . 7 p} (\x] is the ceiling function 

of x). The latter statement is proved in Proposition 13.31 below. It is possible to visualize these Casimirs and 
integrals of motion in the form of the following triangle 

H ^m.-\- 2 ^[^l.o H \™],i H \ g f]' 2 ^ 3 ' 36 ^ 



Hn,-p • ■ ■ ■ • • • Hnp 

If N and p are relatively prime, Casimir functions Hn,± p are proportional to those introduced in p. 2511 

N 
i=l 

Also, one can note that -Hjv.p, if/v.-p and the second order Casimir functions «SVi+p<Si+p.i, i G {1, ■ ■ ■ ,N}, 
are not independent. 

Proposition 3.3. The coefficients Hk(j\±j, k(j) = \jN/p\ , j € of expansion $3. 33\) are Casimir 

functions. 

Proof. Functions Hk(j),±j are the coefficients of the terms w l 3 = e (±jzi) of the expansion of Tr I L rot (z\) J 
in (|3.32p . To prove the statement of the proposition we will show that terms w 1 3 are not present in the expansion 
of |t- (Z rot (zi)) ,Z rot (z2)| as a series in wi. Using representation (|3~26| of linear brackets (|3~23| . (|3~24]) 
and the following form of r-matrix 



A — ' L n J 

m.n 

where 

(7T777 \ 
-»rj ~ cot ( 7r ?) 7 n = 0, < 771 < N, 



sin(7r2) \2 N 



7T / Z TIZ 

: e [7; ~ "TF ) > < n < N, 0^m<N, 



we obtain 
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Z2) 



T mn , (L™{zx) 



/ —1 

E £ [ m ] (Fi - ?2) £ ( £rot (*)) T ™ (£ rot (^i) 



k—i 



>T_ 



This gives 



fe-1 



(3.37) 



Substituting the explicit expressions of Sin- Algebra basis elements T m „ and 93 [ m l (z) into formula (|3.37p 

we get 



where 



JV-1 N 



| Tr (z- ( i l} ) fc , }=^Ecot e fe) (z- ( 

^ - 1 m=l i=l 

x ((Z™^))^ 

+7rifciV (Z™^)) E f ( £r0t (^) 



Z 2 j ] X 



/ iim \ / Jim 
6 " e "IT 



i=l 



fe-1 



(Sat -S in )~2TrikN- 



w 1 



Wx - w 2 



K,L Iot (z 2 ) , (3.38) 



K ilh = e I - 



{{h ~ ii) mod (zi - z 2 ) 



i('2 
« ! 1 



(il-Jl) mod W 



N 



U°\zx) 



(Z™\zxj) 



fc-1 



fe-1 



By the properties of T n 



II T « 



N 



5 (^1 + E n ' ~ ^ J 



Therefore, each element of matrix _K~ is a Laurent polynomial in wx, which can be seen from expression (|3.33[) . 
More precisely, 



^=E- E 

mi,m m fc _ 1 ,n fc _ 1 
/ fe-1 \ fe-1 

x<5 iix + E 711 ~ 3i I H c ' (™>i,ni) 



('l-Jl) ™° d W ('1 -Jl) mod JV »» 



{=1 



2=1 



and the degree of wx is the integer 



(ix - ji) mod TV _ y4 rt; 



(ii — ji) mod AT xTV + (ji — ii) mod N 



N 



N 



N 



N 



= — [x + 1 — 6^), xeZ. 

To derive the maximum x m „„ and minimum x m -, n of x we use the condition 



(3.39) 
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k-1 

1=1 



= \xN + (ji - ii) mod N\ < (k - l)p. 



Consequently, 



(k — l)p — (J i — ii) mod N 
N 



(k — l)p + (ji — ii) mod N 



N 



In the case when k = \jN/p\ we have [_(k — l)p/N\ = j — 1 and hence 



Moreover, x m i n = 1 — j for the diagonal elements of matrix K. According to formula (|3.39p we obtain the 
maximum and minimum degrees d max and d mm of W\ in matrix K, respectively, 



Hence, each element of matrix 



^min 3t ^max j L 

K, L mt (z2) is a Laurent polynomial in w\ with maximum and minimum degrees 



and dmin. Due to the fact that 



polynomials w{ K,L rot (z2) 



K\ W1=W2 = (L rot (z 2 ) 
have root w\ — W2 and 



fc-i 



-2ttLV- 



Wi 



Wi - w 2 



K,U ot {z 2 ) 



(the last term in (|3.38p ) is a Laurent polynomial in w\ with degrees from the interval [1 — j,j — 1]. Also, the 

other two terms in f|3 . 38[) contain only the diagonal elements of matrix I U ot {z\ ) j and thus they are Laurent 

polynomials in u>\ with degrees from the same interval [1 — j,j — 1]. Therefore, terms w^ 3 are not present in 
the expansion (|3.38[) . □ 

Furthermore, it is convenient to treat Hkj as functions of the following set of variables 



S=L, ie{l,...,N}, £> = 0; S jk , j^k, j,ke{l,...,N}j, 

where there are N — 1 independent coordinates ij and N — 1 independent coordinates Sjj , wich are connected 
through nondegenerate transformation (|3.29l) . 

Proposition 3.4. If N and p are relatively prime, then the coefficients Hkj, k > 0, \_jN/p\ ^ k ^ N, \j\ < p, 
the second order Casimir functions Si^+ p Si+ p ^, 1 ^ i ^ N, and H^ p are functionally independent. 

Proof. Coefficients Hkj are the fc'th order homogeneous polynomials in variables S mn . To prove the indepen- 
dence we are going to consider the terms with the maximum degree of the variables {li, 1 ^ i ^ N}. These 

terms in turn contain ones with the maximum degree of the variables ^S^iip, 1 ^ i ^ A |. Let us denote these 

terms by H' k j. Then the independence of Hkj follows from the independence of H' k j. Moreover, the statement 
of the proposition follows from the independence of H' k j, k > 0, [jN/p\ ^ k ^ N, \j\ < p, and the second 
order Casimir functions (the second order Casimir functions are monomials and remain the same after taking 
terms with the maximum degree of any set of variables). After taking the leading terms in Casimirs we obtain 

N 

H X0 « X^ 4 ' 



N k(j)-l 

e n s 

i—l \ rn — 1 



i±(m--l)p,i±mp j $izk(k(j) — l)p,i 5 



m = 



Nj 



0<j<p. 



(3.40) 



By the properties of Casimirs as the coefficients in (|3.32l) . H' k j are homogeneous polynomials in li's and the 
summands of (|3.40[) . so it is convenient to treat these summands as new variables 
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x„ 



i±(k(j)-l)p,i 



'kGO-1 

n * 

m— 1 



i±(m-l)p,i±mp I ^i±(k(j) — l)p,i > 



The expression 



<Si±(fc(j)-l)p,i — ^[^JiM — ^iiiVj+W) mod p,i 



1 < i < AT, < j < p. 



^z=F(./Vj) mod p,'. 



implies the independence of X in the case when N and p are relatively prime and < j < p. 
Therefore, we introduce the set of variables 



S 



1 = {k, 



1 < i < N: 



Xi, i±j , (1 < i SC N) A (0 < j < p); S M±P , 1 < i SC iv} 



and expand the differentials of the functions H' k -, 1 < k ^ N, \j\ < P, H' Np , and the second order Casimirs 
in the basis of the differentials of S 1 . We treat the differentials of li, 1 ^ i ^ N, as independent due to the 
presence of H' w . which determines the connection between them. That leads to the Jacobian matrix of the map 
from S 1 to Hamiltonians and Casimirs. The independence of functions under consideration follows directly from 
the fact that the Jacobian matrix has maximum rank. After appropriate ordering of variables, Hamiltonians, 
and Casimirs one can write the Jacobian matrix in a block lower triangular form 



ff'(0) 



H'(j) 
H'(-j) 



TTl 

n Np 



\ 



■ ■ In 




x(j) 


x(-j) 




SpN 


^ $i.i-\-p 


V 











































A{j) 




















m 










































b 






















D 



(3.41) 



Since we care about the rank of the Jacobian matrix, without loss of generality we can neglect numerical 
common factors in each H' k ■ and use the ordering 



H'(Q) 



H'(j) 




< |j| < p; 



\ Sn,n+ p Sn+ p ,n 



V H N0 

Also, in p.4ip we denote 

X(j) = {Xi_(jvj) mod P ,i, i = l + mp, < m < N - k(j)} , < j < p 



X(-j) = {X iti -( N j) mod P , i = l + mp, < m < N - k(j)} , < j < p, 
and sort out variables in amount equal to the number of elements in H'(±j) so that A(j) are square matrices, 

{S iti+P } = {S iii+P , 1 < i < AT}. 

Due to the fact that every diagonal block of the matrix (|3.4ip is square, we can calculate the determinant 
of p. 4111 . The first diagonal block is a square Vandermonde matrix 



V = 



1 I 



1 \ 

In 



N-l 



L N / 
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with the well-known determinant det V = JJ (lj — U). 

The blocks A(j) are more complicated. Every element A mn (j) of A(j) is a complete homogeneous symmetric 
polynomial (see [22]) h m -i(n,j) in {Z 1+ ( n _ 1)p+rp , < r < - l}, i.e., 

+00 fc(j)-i 
A mn {j) = h m ^i{n,j), ^ h m (n,j)t m = ]Q (l - Ui + ( n _i) p+rp ) . 

m— r — 

The determinants of these matrices are 

det AO') = J] ^ 

il>*2 

ii = {l + mp, < m < 2 (fc(j) - 1)} , i 2 = {l + mp, < m < fc(j)} , 

(|E.3j) (see|E|. The last two diagonal blocks in ()3.4ip are simple: 

dH' N dH Np YJ o 

= ^r— - = — ^— 51 OC I I 

= diag l^i+j,,!, 5 , 2+ p ,2, • • • , Spjv j • 
Thus, the determinant of the matrix (|3.41[) is the following product: 

p-i 

bdetV dct D]J (detA(j)) 2 
and the Jacobian matrix has the maximum rank. 

□ 

Now we are to prove the main statement of this subsection. 

Proposition 3.5. If N andp are relatively prime, then the systems under consideration are completely integrable 
in the Liouville sense. 

Proof. Poisson algebra of the systems can be separated into two disjoint subalgebras. As it was stated earlier in 
Subsection 13.3.11 the elements of the first subalgebra are the functions of the elements of the second one on the 
generic symplectic submanifold. The submanifold corresponding to the elements of the second subalgebra has 
the dimension (2p + 1) N — 1. There is the following condition (see[Dj for the dimension R of the symplectic 
leaf of this submanifold: 

R > 2p(N- 1). 

Propositions 13.31 and 13.41 give us 2p + N — 1 independent Casimir functions and, consequently, we have the 
equality 

R=2p(N-l). 

According to the Liouville theorem about integrable systems it is necessary to have R/2 functionally inde- 
pendent Hamiltonians in involution for the complete integrability. From Propositions 13.21 13.41 and (|3.36p we 
have p(N — 1) = R/2 independent Hamiltonians in involution. 

□ 



4 Conclusion 

We have proposed a procedure giving a limit relation between the elliptic SL(N, C) top and the Toda chains. 
This procedure is similar to the Inozemtsev limit and is a combination of the shift of the spectral parameter, 
the scalings of coordinates and the trigonometric limit. 

Also, in Subsection 13.31 we have shown, that the generalization (|3.2ip of the above procedure provides a 
new class of integrable systems in the case when N and p > 1 are relatively prime. The open problem is to 
understand whether the limiting systems are integrable in general case when p < N/2. Some statements, such 
as Propositions 13.11 13.21 13.31 are still valid in general case, but the whole set of independent Casimir functions 
and Hamiltonians is not clear. 
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A Sin-Algebra 

We will use the following notation to simplify formulae: 



S(n) = 



1, n = mod N, 
0, n ^ mod N, 



e(z) = exp (27riz) . 

Elements T mn of the Sin- Algebra basis in si (N, C) can be written in the form 
/ mn \ I iTfi \ — 

Crm») y =e(— ]e( — )8(j-i-n), m^Q or n + 0, m, n € {0 . . .N - 1}. 

For m,ngZ, (m ^ mod N) or (n ^ mod N) , the quasi-periodic condition can be introduced 

/ mn — (m mod N)(n mod N)\ 

' mn — (> I 2JV / mod N ,n mod JVi 



(m mod N)(n mod A) — rrm 
2A 



mod N.n mod jVi 



where e ((mn — (m mod A)(n mod N)) J (27V)) = ±1. 
The commutator relations in this basis are 

r it i 

[T mn ,T k i] = 2isin (kn - ml)^ T m+k>n+ i. (A.l) 

One can establish the following relations between the coordinates in the standard {SV,} and Sin- Algebra 
{s mn } bases 

$ij — ^ ^ Srnn (^mn) ij ■> ^mn — ~T7 ^ ^ &%j (^— m,n)^- ■ (-^ - ^) 



N 



B Degenerations of elliptic functions 



Definitions and properties of elliptic functions are borrowed mainly from |11) and |23j . The main object is the 
theta function with characteristics defined via 



(z,T) = ^0-+") 2 e(( J+a )( z + 6)). 



where q — e (r) = exp (27rir). 

We will also need the Eisenstein functions 

M M 

e k (z)= lim V (z + n)~ k , E k (z) = lim V e k (z + nr). (B.l) 

M— >-|-oo * — ' M— >+oo » — ' 

n=-M n=— M 

To determine the limits of Lax matrices we will use the series expansions of the following functions 



d(z) 



1/2 
Ll/2 



,r) = ^S j+ ^e((j + l)(z + l)), (B.2 



,/6- 
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(u,z) 
m + tit 



tf(u + z)$'{0) 
■d(u)d(z) ' 



(B.3) 
(B.4) 



The functions satisfy the following well-known identities: 

(j>(u,z)4>(-u,z) = E 2 {z) - E 2 (u), 

d u (j)(u, z) = 4>(u, z){E\(u + z)- Ei(u)), 



(B.5) 



parity 



and quasi-periodicity 





E k (-z 


) = (-l) k E k (z), 






n- 


■z) = -d{z), 






z) = (j) 


z,u) = -<t>{-u,-z), 




£1(2 + 1) 


= E 1 (z), 


E 1 {z + r)=E 1 {z)- 


27ri, 


E 2 (z + l) 


= E 2 (z), 


E 2 (z + t)=E 2 {z), 




§(z + l) = 




i9(z + t) = -q~^e(- 


■z)*{z), 



(B.6) 



4>{u + 1, z) = <p(u, z), <fi(u + t, z) = e(—z)4>(u, z). 
We will examine degenerations of elliptic functions ()B.4|) in the following limit: 

~ T 

z = z + — , Im(j) — > +00. 

Using definition (|B.3[) one can reduce the expansion of ip [™] (z) to the expansion of theta functions. Con- 
sidering the main non- vanishing terms we have 



/to n 



rt / t to n 



= < 



2qs sin^TT^^) +0 (gs) , n = 0, 

(~2^) ^) ' 0<n<iV ' 

^~*e Q ^ +0 (gw-l) , 0a<| 

2gssin(7r(z-^)) +o(g*) , n = y, 

^( 1 -|)) +0 ( ?§ ^)' 



jgrs 2jve 



iqs 2«e 



<n< N, 



wich gives 



to + nT _ r 
iV - ;Z ' 2 



- = < 



-7re f — 1 sin -1 (V— 1 + o (1) 
2itiq%e (J^j +o(gs) , 
4tt<73 sin (tt ^)) 
— 27ri<7^e (z) + o (q 



TO 1_ 

e [ h -z 

2iV 2 



o g 



71 = 0, 



N 

0<n<—, 



N 



(B.7) 



2 ' 



iV 



< n < N, 
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(in 


-i) 


. n - 







-e(— jam- (tt- j + o (1) , 



r. • ^_ / Tl HZ 1 / _n_N 

27Tig2N e ( — - — ) + (( ? 2ivj , 



, i / m \ . ( /_ m 
47T(? 4 e sin i z 



o 



-27rig 2n e 



N-n, 



N 



Z + o I q aw 



n = 0, 



< n < —, 



N 



2 ' 



N 

— <n < N. 
2 



(B.8) 



To evaluate the limit of / [™] we expand E\ (x — ar) as a series in q. From the definition (|B.1|) we get 



A I 



Ei (x — ar) = lim ej (x + (n — cr) r) = £j (x — err) 



-M 



M 



lim > (ei (x + (n — cr) t) + ej (x — (n + cr) r)) 



M->+oo 



Using the explicit formula for £\{x) from [TT] 



-1 — 2e(x) + o (e(x)) , Irn(x) —> +oo, 



f s if \ A%) + 1 

Sl(x) = 7TCOt (7TX) = 7T1— = 7T1 \ , , n / \ , f I \\ T I \ 

w v ; e(x) - 1 \ 1 + 2c(x) + o(e(x)) , Im(x) +oo, 



one can write down the leading term of £i(x — or) as follows: 

7rcot(7rx) + o (1) , 



Ei (x — ar) = < 



Tri + 27ric/ T e(-x) + o (c/ CT ) , 
7ri + 1-K\q^ (e(— x) — e(x)) + o 
Tri - 2niq 1 - a e (x) + o , 



and using (|B.6[) generalize it to cr 6 M: 

27Ti [crj + 7T cot(7rx) + o (1) . 



cr = 0, 

< a < -, 
1 Z 

f 



{cr} = 0, 



Ei(x - err) = < 



0<{cr}<-, 



2rri [crj + rri + 27ric/^>e(-x) + o (g^>) , 

27ri [crj + 7ri + 27ric/5 (e(-x) - e(x)) + o fqi} , {cr} = -, 

1 

27ri [crj + Tri - 27rig 1 -{ <J >e (x) + o (q 1 '^) , - < {cr} < 1, 



where {cr} is the fractional part of cr. 

To expand d u 4>(u, z)\u=u-<tt m the limit Im(r) — > +oo it is convenient to use formula (|B.5|) . Assuming 
z = z + r/2 we have to consider the following cases depending on the value of cr: 



1. cr = 



Ei (u + I+ T - 



o(l), 



sin 7rw 

-El (u) = 7T COt 7TU + O (1) , 

-7ii + o(l) , 



d u (j)(u, z)\ u= u-ar = SU1~ 2 7TU + O (1) 
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2. < a < | 



Ei I u + z - 



(f>(u — err, z + - J = 2n\qe(— u) + o (q) , 
Ei(u — or) = 7ri + 27rige(— u) + o (1) , 

— 7ri — 27ri<7^~ CT e (u + z) + o (^q 

d u <f>(u, z)\ u= u- aT = 47r 2 ge (-u) + o (g) . 



--air 



3. a = i 



4>(u,z) — Anq 2 sin (tt (it + z)) e ^— — u + -z^ + o ^g 2 ^ , 

£?i (u + 2) = 7r cot (n (u + z)) + o (1) , 
1 



Ex ( 5 - -r ) = Tri + o (1) , 



4. i < cr < 1 



4?r' 



! g2 e (-M) +o (cp) . 



U — (TT, Z + — 



El 



-27rig2 e (F) + o fc^) , 
i?i (u — err) = 7ri — 27ric/ 1_cr e (w) , 
(u + z + Q - crj r) = Tri + 27ric/ ff ~3e (— S — 5) + o (V~^) , 

4 



z)| m=s _ ctt = 4tt 2 c/2 
4TT 2 qe (-u) +o(q) , 



e — u — z + cr — 



4rr 2 g4 (e (— it) — e (it + z)) , a 



1 3 

— < <7 < — , 

2 3 4 



— e (w + (1 — it) r) 



e(2) 



4' 



-4tt 2 c/2 <T e(w + z) 



- < cr < 1. 
4 



Summarizing all the special cases we obtain 

— 7T 2 sin -2 TTU + o (1) , 



d u (j)(u, z)| u =a_ CTT = < 



47r 2 c/e (-2) + o (q) , 
A-K 2 qi [e (— u) — e (u + z)] + o (q 
-47r 2 c/l _cr e (2 + z) + o (qi~ a ^J , 



(7 = 0, 

3 

0< Cr<I> 

3 

- < cr < 1. 
4 



Finally, from (|B.4[) we get 

-tt 2 sin" 2 (vr^) +o(l), 

^ (]?) 6 9 ^ + ° ( ? ^) ' 

(^)- e (-^ + F )] e (-I ? ) gl+0 ( 9§ )' 
-47r 2 e (-^ + zj e f-^ J c/^ 1 -^-) + ^K^w) 



"TO" 










-5) 


. n - 







4rr 2 



71 = 0, 



3iV 

< n < — — , 



3iV 



3iV 

— — < n < N. 

4 
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C Scaling inequality 



We are going to prove the inequality Oj + aj — oti+j ^ 0, Vi,j € Z, for the function 



O', = < 



r fc 

2JV' 
P 



0^k<p^—, 



1 



N-p<k<N, 



12 2iV 
where k = i mod AT. 

Taking into account that asj+jy = ctj, it is enough to consider the case i, j G [—AT/2; N/2], where we have 



2N 



, N < P, 



2iV' f<l*l< T 



It is easy to check that there are the following two cases: 
1. If (\i\>p)V (\j\>p), then 

P 



2. H(|i|<p)A(|i|<p),then 



. i+j < I* 

Q!i+j ' — — ' 



\j\ 



Ot; + ft;. 

2N 2N 2N 3 



D Dimension of the symplectic leaf 

The phase space of the system from Subsection 13.31 is equipped with the following Poisson structure: 



{Su,Sjk} = N(Sji5ik — SikSij), 
{Sij,Ski} — N(SkjSu — SuSkj), 

(0 < (J - i) mod N < p) A (0 < (I - k) mod N ^ p) A (0 < (j + I - i - k) mod N < p), 
{Sij,Ski} — N(SkjSu — SuSkj), 

(N - p < (j ~ mod N <N)A(N-p^(l-k) mod TV < AT) A 
A (AT - p < (j + Z - i - fc) mod iV < TV). 

We are interested in the specific subalgebra 

S2 = {Sij, (0 < {j - i) mod iV < p) or (AT - p < (j - i) mod N < AT)| 
and the restriction of the Poisson tensor 7rwH H ) (<S 2 ) on it 



(5 2 ) , G (5 2 )} = 7r<«>< H > (5 2 ) a (H) G, d { 



d 

OSi 



The dimension of the symplectic leaf of the Poisson submanifold <S> 2 is the rank R of 7rW)w (<S> 2 ). Our aim 
is to minorize the rank R. Note that 7rfa)'*') (<S 2 ) can be represented as a square (2p+ 1)N x (2p + l)AT-matrix. 
To write down this matrix in a block triangular form we use the following ordering: 



Y = {Y(k), k = ±p, . . . ,±1,0} , Y(k) = {Si, t 



+k, i = 1, 



,N 
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Thus, we have 



where 



Y(p) Y(-p) 



Y(p) 
Y(-p) 

Y(p-1) 
Y(l-p) 



Y(l) 
Y(-l) 

Y(0) 



\(P+(p) P-(p)) 



Y(p-l)Y(l-p) 




p+(p - i) o 

P_(p-1) 



Y(l) Y(-l) 


P + (l) 
P_(l) 



Y(0) 

P+(0) 
P_(0) 



(D.l) 



(P + (k))ij — js^i+p^fc, S'jj+fcj — N (^Sj^ + p5 it j + k — Si^ + p5j t i +p -kj , k G {0, . . . ,p} , 

(p_(fe))y = j,s M+fc _p, •S'jj-fcj = n (Sjj ./),., /. - SM-p^+fc-p) , fc e {o, . . . ,p} . 

It can be easily seen that the rank of each matrix P±(k) is equal to N — 1. Due to the block triangular 
form of the matrix (|D. 1|) we obtain the required condition 



R > 2p(N- 1). 



(D.2) 



E Addition to Proposition 13. 4L det A(j) 

Here we are to consider square matrices A(j) with the following structure: every element A mn (j) of A(j) is a 
complete homogeneous symmetric polynomial (see |22| ) /i m _i(n,j) in variables {Zi_|_(„_i)p4. r p, O^r^ fc(j) - 1}, 



i.e. 



+00 



k(j')-l 



4m(j) = h m -i(n,j), ^ h m{n,j)t m = J| (l - ^l + (n-l)p+rp) 1 , 



(E.l) 



m— 



r=0 



where we use the notation k(j) — \jN/p\ from Proposition! 

The determinants of the matrices under consideration are homogeneous polynomials of order 
(k(j)(k(j) — 1)) /2. In order to compute these determinants we will find the appropriate number of roots and 
calculate their values at one particular point. 

As one can see from (|E.1[) columns in matrix A(j) depend on the following set of variables: 



Column Number 

1 l\ /i+p 

2 ^i+p 



Set of Variables 

h+(k(j)-i) P h+k(j) P 



(E.2) 



k(j) ?i + ( fe (j)_i) p h+k{j) P ■ ■ ■ h+2(k(j)-i) P 

It follows from (|E.2|) that two adjacent columns will coincide if we put the first variable, e.g., 1%, in one 



column equal to the last variable, e.g., li+Mfipi in the subsequent column. Thus, we have \jN/p\ — 1 roots. 
It turns out that the result can be generalized to any c adjacent columns, where 2 ^ c ^ k(j) — \jN/p]. 
In other words, if the first variable in the first column is equal to the last variable in the last column, then 
the determinant is equal to zero. Indeed, without loss of generality we can consider the first c columns. By 
assumption, 



i+(fe(i) 



c-2)p — h- 



We will prove that in this case the rows of A(j) are linearly dependent. If we multiply the m'th row of 
matrix by the elementary symmetric polynomial (^^) k ^^ m £k(j)-m(j) m variables 
{h +rp , ^ r < 2k(j) - 2, r ? k(j) + c - 2} 
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2k(j)-2 

m^O r=0,r^k(j)+c-2 

then after summation we get the zero row due to the fact that the elements of the row are coefficients in front 
of t fc 0') _1 in the following polynomial with the maximum degree k(j) — 2: 



m^O mi^O 

(llr=0 (! - 1 il+rp)) (n'=fcfeln-l,r^fe0)+c-2 ( X ~ f *l+rp)) . " < C > 

= < 

k (to (i - 1 h +rP )) (n'Sfe+n-i (! - 4 ^)) . n > c - 

Thus, we have (k(j)(k(j) — 1)) /2 roots. Now we calculate the determinant at the following point 

/ m + n — 3 

Hi — t, ii 2 — U, ^-ran I 

\ m — 1 

ii = {l + mp, fc(j) < m < 2 (fc(j) - 1)} , i 2 = {l + mp, < m < fc(j)} . 
Subtracting columns according to the well-known formula, i.e., 

'N + 1\ _ /A\ / TV 
m+1/ \m/ \m+l / 

one can obtain a lower unit triangular matrix, so at the point selected above we have 

dct A = 1. 

Finally, we get 

dot - n {h.-hi), 

H>»2 

^l = {l + mp, m < 2 (fc(j) - 1)} , i 2 = {l + mp, < m < fc(j)} . (E.3) 
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